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ABSTRACT 

We present full evolutionary calculations appropriate for the study of hot hydrogen-deficent DO 
white dwarfs, PC 1159 stars, and DB white dwarfs. White dwarf sequences are computed for a wide 
range of stellar masses and helium envelopes on the basis of a complete treatment of the evolutionary 
history of progenitors stars, including the core hydrogen and helium burning phases, the thermally- 
pulsing AGB phase, and the born-again episode that is responsible for the hydrogen deficiency. We 
also provide colors and magnitudes for the new sequences for Toff < 40 000 K, where the NLTE 
effects are not dominant. These new calculations provide an homogeneous set of evolutionary tracks 
appropriate for mass and age determinations for both PG 1159 stars and DO white dwarfs. The 
calculations are extended down to an effective temperature of 7000 K. We applied these new tracks to 
redetermine stellar masses and ages of all known DO white dwarfs with spectroscopically-determined 
effective temperatures and gravities, and compare them with previous results. We also compare for 
the first time consistent mass determinations for both DO and PG 1159 stars, and find a considerably 
higher mean mass for the DO white dwarfs. We discuss as well the chemical profile expected in the 
envelope of variable DB white dwarfs from the consideration of the evolutionary history of progenitor 
stars. Finally, we present tentative evidence for a different evolutionary channel, other than that 
involving the PG 1159 stars, for the formation of hot, hydrogen-deficient white dwarfs. 
Subject headings: stars: evolution — stars: abundances — stars: AGB and post-AGB — stars: 
interiors — stars: variables: other (GW Virginis) — white dwarfs 



1. INTRODUGTION 

White dwarf stars are the end-product of the evolution 
of the vast majority of stars. Indeed, more than 97% of 
all stars in our Galaxy are expected to end their lives 
as white dwarfs. For this reason, the present population 
of white dwarfs contains valuable information about the 
evolution of individual stars, the previous history of the 
Galatic populations (Isern et al. 1998; Torres et al. 2002) 
and the rate of star formation (Dfaz-Pinto et al. 1994). 
In addition, white dwarfs have potential applications as 
reliable cosmic clocks to infer the age of a wide variety of 
stellar populations (Winget et al. 1987; Garcfa-Berro et 
al. 1988; Hansen et al. 2007), and to place constraints 
on elementary particles (Corsico et al. 2001; Isern et al. 
2008). Recent reviews on the properties and evolution 
of white dwarfs and of their applications are those of 
Winget & Kepler (2008) and Fontaine & Brassard (2008). 

Traditionally, white dwarfs have been classified into 
two distinct families according to the main constituent 
of their surface: those with a hydrogen (H)-dominated 
atmosphere — the DA spectral type — and those with 
a helium (He)-rich surface composition — the non-DA 
white dwarfs. The latter, that comprise about ^ 20% 
of known white dwarfs, are usually divided into three 
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subclasses: the DO spectral type (with effective temper- 
atures 45 000 K < Toff < 200 000 K) that shows relatively 
strong lines of singly ionized He (Heii), the DB type 
(11000 K < Tefi < 30 000 K), with strong neutral He 
(Hei) lines, and the DC, DQ, and DZ types (Tofi < 11 000 
K) showing traces of carbon and metals in their spec- 
tra. As a DO white dwarf evolves, the Hen recombines 
to form Hei, ultimately transforming into a DB white 
dwarf. The transition from DO to the cooler DB stage 
is interrupted by the non-DA gap (that occurs at 30 000 
K < Toff < 45 000 K) where few objects with H-deficient 
atmospheres have been observed (Eisenstein et al. 2006). 
To this list, we have to add the discovery of a new white 
dwarf spectral type with carbon-dominated atmospheres, 
the "hot DQ" white dwarfs, with Toff - 20 000 K (Dufour 
et al. 2007, 2008). 

DO white dwarfs are usually thought to be the progeny 
of PG 1159 stars (Dreizler & Werner 1996; Unglaub 
& Bues 2000; Althaus et al. 2005), hot stars with H- 
deficient and He-, C- and 0-rich surface layers (Werner 
& Herwig 2006). A large fraction of PG 1159 stars is 
believed to be the result of a born-again episode, i.e., a 
very late thermal pulse (VLTP) experienced by a white 
dwarf during its early cooling phase (Fujimoto 1977; 
Schonberner 1979; Iben et al. 1983; Althaus et al. 2005). 
During the VLTP, most of the H content of the remnant 
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is violently burned (Herwig at al. 1999; Miller Berto- 
lami et al. 2006). As a result, the remnant is forced 
to evolve rapidly back to the AGB and finally into the 
central star of a planetary nebula as a hot H-deficient 
object and with a surface composition resembling that 
of the intershell region chemistry of AGB star models — 
typical mass abundances are AThc — 0.33, Xq — 0.5 and 
Xo — 0.17, though notable variations are found from star 
to star (Dreizler & Heber 1998; Werner 2001). Eventu- 
ally, gravitational settling acting during the early stages 
of white dwarf evolution causes He to float and heavier 
elements to sink down, giving rise to an He-dominated 
surface, and turning the PC 1159 star into a DO white 
dwarf (Unglaub & Bues 2000). During this stage, the 
evolution of the star is dictated essentially by neutrino 
losses and the release of gravothermal energy (O'Brien & 
Kawaler 2000; Althaus et al. 2005). Because these white 
dwarfs are the hottest ones, the evolution during the DO 
stage proceeds very fast, with typical evolutionary time 
scales ~1 Myr. 

Several works have been devoted to the identification 
and spectral analysis of H-deficient objects like PG 1159 
and DO white dwarfs — see Dreizler & Werner (1996), 
Dreizler et al. (1997) and references therein for earlier 
studies of spectroscopically confirmed DO white dwarfs. 
Particularly relevant is the star KPD0005-I-5106, the 
hottest known DO white dwarf with an effective tem- 
perature of about 200 000 K (Werner et al. 2008a), 
much higher than previously thought (Dreizler & Werner 
1996). With the advent of large surveys, like the Sloan 
Digital Sky Survey (York et al. 2000) the number 
of spectroscopically-identified H-deficient stars has in- 
creased considerably. For instance, Krzesihski et al. 
(2004) reported 15 spectroscopically-identified DO stars 
— 13 of which were new discoveries — from the SDSS 
Data Release 1 (Abazajian et al. 2003), which contains 
«2 500 white dwarfs (Kleinman et al. 2004). Also, Eisen- 
stein et al. (2006) reported 31 DO white dwarfs and 
10 PG 1159 stars, from their catalog based on SDSS 
DR4 (Adelman-McCarthy et al. 2006). Hiigelmeyer 
et al. (2005) performed spectral analyses of 16 hot H- 
deficient stars from DRl, DR2 and DR3 catalogs , where 
9 were classified as DO and 5 as PG 1159 stars. More 
recently, Hiigelmeyer et al. (2006) extended this work 
to 13 DO and 4 PG 1159 stars — see also Werner et al. 
(2004). These authors derived effective temperatures, 
surface gravities and stellar masses for their sample, and 
for that of Hiigelmeyer et al. (2005), and discussed the 
evolutionary connection between DO white dwarfs and 
PG 1159 stars. 

To infer the mass and age of observed DO white dwarfs, 
evolutionary calculations for these stars are required. 
However, most of the existing calculations that are em- 
ployed in the determination of evolutionary parameters 
of non-DA white dwarf stars (Wood 1995; Benvenuto & 
Althaus 1999) are not entirely consistent with the evolu- 
tionary history that leads to the formation of H-deficient 
stars. Even though significant deviations are not ex- 
pected at low effective temperatures, this lack of con- 
sistency becomes relevant for the high effective tempera- 
tures characteristic of hot DO white dwarfs. In addition, 
the lack of evolutionary calculations covering the whole 
stage from the domain of luminous PG 1159 to the hot 
white dwarf stage prevents from a consistent mass de- 



TABLE 1 

Initial and final stellar, mass (in solar, units), and the 

TOTAL MASS OF He CONTENT LEFT IN THE WHITE DWARF FOR THE 
EV0LUTI0NAR,Y SEQUENCES CONSIDERED IN THIS WORK. 





-'^^ZAMS 




0.515 


1.00 


0.0219 


0.530 


1.00 


0.0090 


0.542 


1.00 


0.0072 


0.565 


2.20 


0.0067 


0.584 


2.50 


0.0060 


0.609 


3.05 


0.0058 


0.664 


3.50 


0.0036 


0.741 


3.75 


0.0019 


0.870 


5.50 


0.0009 



termination for the DOs and PG 1159 stars (Dreizler & 
Werner 1996). To our knowledge, the only exception is 
the study of Althaus et al. (2005), who computed the for- 
mation and evolution of a 0.588 Mq white dwarf remnant 
from the ZAMS through the born-again and PG 1159 
stages to the white dwarf regime. This paper is precisely 
intended to fill this gap, by computing new evolution- 
ary tracks and mass-radius relations appropriate for both 
non-DA white dwarfs and PG 1159 stars. Sequences are 
derived from a full and self-consistent treatment of the 
complete evolutionary history of progenitor stars with 
different masses evolved through the born-again episode. 
In this way, the white dwarf evolutionary calculations 
presented here are not affected by inconsistencies aris- 
ing from artificial procedures to generate starting white 
dwarf configurations. In addition, the computation of 
the evolutionary history of progenitor stars gives us the 
amount of He left in the white dwarf as well as the chemi- 
cal profiles expected not only for the carbon-oxygen core, 
but also for the partially degenerate regions above the 
core, of relevance for the white dwarf cooling phase. 

The calculations presented here constitute an homoge- 
neous set of evolutionary tracks aimed to study both PG 
1159 and DO white dwarf stars, as well as the DB white 
dwarfs. The paper is organized as follows. In Sect. [2] we 
describe the main physical inputs to the models. We also 
describe the main characteristics of the initial models and 
some details of the evolutionary computations. In Sect. 
Owe present the results for DO white dwarf evolutionary 
sequences. We apply the sequences to redetermine stellar 
masses and ages of observed DO white dwarfs, and we 
derive mass-radius relations. In addition, we show the 
results for the case of DO white dwarf sequences with 
very thin He-rich envelopes. We close the paper in Sect. 
[3] by discussing and summarizing the implications of our 
results, particularly regarding the origin of the observed 
DO white dwarfs. 

2. COMPUTATIONAL DETAILS 

2.1. Input physics 

The evolutionary calculations presented in this work 
were done with the LPCODE stellar evolutionary code we 
employed in our previous studies on the formation of 
the H-deficient stars through late He flashes, like PG 
1159 and extreme horizontal branch stars (Althaus et al. 
2005; Miller Bertolami & Althaus 2006; Miller Berto- 
lami et al. 2008). More recently, the code was used to 
compute the formation of hot DQ white dwarfs (Althaus 
et al. 2009b) as well as the formation and evolution of 
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He-core white dwarfs with high metaUicity progenitors 
(Althaus et al. 2009c). Details of LPCODE can be found 
in these works. In particular, the code considers a si- 
multaneous treatment of non-instantaneous mixing and 
burning of elements, which is of primary importance for 
the calculation of chemical abundance changes during the 
short-lived evolutionary stages characteristic of unstable 
burning episodes, like the born-again stage, from which 
our starting H-deficient white dwarf configurations are 
derived, see 2.1. Nuclear reaction rates are from Caugh- 
lan & Fowler (1988) and Angulo et al. (1999), and are 
detailed in Althaus et al. (2005). In particular, the 
^^C(a,7)^^0 reaction rate was taken from Angulo et al. 
(1999), which is about twice as large as that of Caughlan 
& Fowler (1988). A moderate, diffusive overshooting in 
the core and in the envelope is allowed during pre-white 
dwarf evolution. 

For the white dwarf regime, we considered the follow- 
ing main physical ingredients. Neutrino emission rates 
for pair, photo and bremsstrahlung processes are those 
of Itoh et al. (1996). For plasma processes we included 
the treatment presented in Haft et al. (1994). Radia- 
tive opacities are those of the OPAL project (Iglesias & 
Rogers 1996), including carbon and oxygen rich compo- 
sition. We adopted the conductive opacities of Cassisi 
et al. (2007). This prescription covers the whole regime 
where electron conduction is relevant. For the high den- 
sity regime, we used the equation of state of Segretain et 
al. (1994), which accounts for all the important contri- 
butions for both the liquid and solid phases — see Al- 
thaus et al. (2007a) and references therein. For the low- 
density regime, we used an updated version of the equa- 
tion of state of Magni & Mazzitelli (1979). The release 
of latent heat upon crystallization was also considered. 
In particular, for our calculations, crystallization sets in 
when the ion coupling constant reaches F = 180, where 
r = (Z5/3)eVaefcBT and fle is the interelectronic dis- 
tance. We considered a latent heat release of fceT per ion, 
which was spread over the range 175 < F < 185. Con- 
vection was treated in the formalism of the mixing length 
theory as given by the ML2 parameterization (Tassoul et 
al. 1990). 

All our white dwarf sequences were computed in a con- 
sistent way with the evolution of the chemical abundance 
distribution caused by element diffusion along the whole 
cooling phase. In particular, we considered gravitational 
settling and chemical diffusion of "'He, ^^C, '^'^C, '^''N and 
'^^O — see Althaus et al. (2003) for details. In our sim- 
ulations, we obtain high amounts of '^'^C and ^^N at the 
VLTP, in agreement with observational inferences in the 
Sakurai object and, in the case of "'"'N, in some PG 1159 
stars (Werner & Herwig 2006). Our treatment of time- 
dependent diffusion is based on the multicomponent gas 
treatment presented in Burgers (1969). Diffusion veloci- 
ties are evaluated at each evolutionary time step. For the 
white dwarf evolution, is worth mentioning that opacity 
is calculated taking into account the diffusion predictions 
for the heavy element composition. 

2.2. Initial models 

The initial models for our white dwarf sequences corre- 
spond to realistic PG 1159 stellar configurations derived 
from the full evolutionary calculations of their progeni- 



tor stars for solar metaUicity (Miller Bertolami & Althaus 
2006). All the sequences were computed from the ZAMS 
through the thermally-pulsing and mass-loss phases on 
the AGB and finally to the born-again stage where the 
remaining H is violently burnt. After the born-again 
episode, the H-deficient, quiescent He-burning remnants 
evolve at constant luminosity to the domain of PG1159 
stars with a surface chemical composition rich in He, car- 
bon and oxygen (Miller Bertolami & Althaus 2006). This 
new generation of PG 1159 evolutionary models has suc- 
ceeded in explaining both the spread in surface chemical 
composition observed in most PG 1159 stars and the lo- 
cation of the GW Vir instability strip in the log Tag— log g 
plane (Corsico et al. 2006), as well as the short born- 
again timescales of V4334 Sgr (Miller Bertolami et al. 
2006). To our knowledge, this is the first time that such 
a detailed treatment of the evolutionary history of pro- 
genitors stars with different initial stellar masses is taken 
into account in the calculation of cooling sequences for 
H-deficient white dwarfs. 

Specifically, we computed nine H-deficient white dwarf 
evolutionary sequences. In Table [T] we list the stellar 
masses of the resulting white dwarfs, together with the 
inital masses of the progenitor stars at the ZAMS. Also 
listed in Table [1] is the mass of He left by previous evo- 
lution. Note that for these sequences, which cover most 
of the observed stellar mass range of H-deficient white 
dwarfs, the expected range of He envelope masses spans 
more than one order of magnitude. Due to He burning 
during the PG 1159 stage, the mass of He left when the 
remnant reaches its cooling track is considerably lower 
than that obtained shortly after the born-again episode 
(0.0036 and 0.0075 M©, respectively, for the 0.664 M© 
sequence). For the 1 M© model, two different sequences 
were computed with two different mass-loss rates during 
the AGB evolution. In this way we obtain a different 
number of thermal pulses and, in the end, two different 
remnant masses of 0.530 and 0.542 M©, respectively. Fi- 
nally, the sequence with 0.515 Af© was specifically com- 
puted to reproduce the observational features of low- 
mass. He-enriched PG 1159 stars, like MCT 0130-1937 
and HS 1517-f 7403 — see Miller Bertolami & Althaus 
(2006) and Althaus et al. (2007b). This model was ex- 
tracted from the full evolution of a 1 M© model star that 
experiences its first thermal pulse as a late thermal pulse 
(LTP) after leaving the AGB. We adopt the mass-loss 
rates of of Blocker (1995) multiplied by a 0.2 factor. Mass 
loss during the departure from the AGB has been arbi- 
trarily set to obtain a final helium shell flash during the 
early white dwarf cooling phase. Our calculations cover 
the hot pre-white dwarf stages (PG 1159 regime) and the 
advanced phases of evolution down to an effective tem- 
perature of 7000 K, therefore covering all the DB cooling 
sequence. 

In Fig. [1] we show the mass abundances of ''He, '^^C 
and ^^O throughout the deep interior of two PG 1159 star 
models corresponding to the sequences with 0.515 and 
0.870 M©. The models shown correspond to evolution- 
ary stages at high luminosities where the chemical profile 
at the bottom of the He-rich envelope is being modified 
by residual He shell-burning. Note the dependence of the 
composition profile on the stellar mass, which emphasizes 
the need for a detailed knowledge of the progenitor his- 
tory for a realistic treatment of white dwarf evolution. 
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TABLE 2 

Main parameters associated to the model with highest Tcff for each sequence. The tabulated abundances (by mass 

fraction) correspond to the surface value. 



M/Mq 


logT.ff [K] 


log{L/LQ) 


log g [cm s ^] 




Xc 


Xo 


0.515 


5.0634 


2.6868 


6.6713 


0.7437 


0.1637 


0.0279 


0.530 


5.1366 


2.9762 


6.6876 


0.3301 


0.3944 


0.1716 


0.542 


5.1650 


3.0546 


6.7412 


0.2805 


0.4064 


0.2127 


0.565 


5.2105 


3.1565 


6.8304 


0.3632 


0.2737 


0.2173 


0.584 


5.2398 


3.2574 


6.8615 


0.3947 


0.3060 


0.1704 


0.609 


5.2819 


3.3533 


6.9518 


0.4990 


0.3465 


0.1033 


0.664 


5.3578 


3.3611 


7.0813 


0.4707 


0.3260 


0.1234 


0.741 


5.4535 


3.8069 


7.2701 


0.4795 


0.3361 


0.1390 


0.870 


5.5829 


4.0961 


7.5679 


0.5433 


0.3012 


0.0938 
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Fig. 1. — Abundance by mass of *He, ^^C and ^''O as a func- 
tion of the outer mass fraction log{l — Af^/M,) for PG 1159 
models corresponding to the sequences of 0.870 (top panel) and 
0.515 Mq (bottom panel). Luminosity and effective temperature 
(log(I//I/0); logTeff) are specified for each model. 



The chemical stratification of our H-deficient pre-white 
dwarf models consists of a carbon-oxygen core, which 
emerges from core He burning in prior stages. The core 
chemical profile is typical of situations in which extra 
mixing episodes beyond the fully convective core during 
the core He burning, like moderate overshooting and/or 
semiconvection, are allowed — see Straniero et al. (2003) 
and also Prada Moroni & Straniero (2007) for the con- 
sequences on white dwarf evolution. The core is sur- 
rounded by a homogeneous and extended envelope rich 
in He, carbon and oxygen, which is the result of prior 
mixing and burning events during the thermally pulsing 
AGB and born-again phases. The surface abundance of 
He of our pre-white dwarf models is within 0.28 and 0.75 
by mass (see Table [2]), which is in agreement with the 
range of observed He abundances in most PG1159 stars 
(Werner & Herwig 2006). These abundances are not only 
determined by the stellar mass, but also by the number of 
thermal pulses during the AGB phase, by the inner pen- 
etration of the convective envelope and by the mass lost 
after the complete burning of protons (Miller Bcrtolami 
& Althaus 2006). Finally, we have removed any H from 
our PG 1159 models, despite that some H is expected to 
be present in PG 1159 stars. This is particularly true if 
some of these stars result from a LTP episode. Here, the 
H content with which the progenitor star departs from 
the AGB for the first time is not burnt during the late 
helium flash but instead diluted to deeper layers. In this 
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Fig. 2. — Upper panel: Hertzsprung-Russell diagram for our evo- 
lutionary sequences. Bottom panel: evolution in the log T^ff — log g 
plane. Filled dots on each track indicate the instant with the high- 
est effective temperature, which we adopt as our time origin. Mass 
values are given in solar units. 

case, the presence of H in PG 1159 stars may be respon- 
sible for the fact that not all PG 1159 stars evolve to 
DO white dwarfs, see Sect. H) By contrast, during the 
VLTP, only traces of H survive the violent H burning, 
see Miller Bertolami et al. (2006). It is expected that 
a large fraction of this remaining H is removed by mass 
loss during the following return to the giant state and 
Sakurai state. 

3. EVOLUTIONARY RESULTS 

3.1. Sequences with thick envelopes 

In Fig. [21 we show the evolution of our H-deficient se- 
quences both in the Hertzsprung-Russell diagram and in 

the log Tcff log 5 plane. The evolutionary sequences 

shown cover from the PG 1159 stage at intermediate 
gravities to the white dwarf domain. Thus, these se- 
quences constitute a homogeneous set of tracks that con- 
sistently cover both evolutionary stages. Some charac- 
teristics of them at the highest effective temperature are 
listed in Table H 
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Fig. 3. — Internal abundance distribution of *He, ^^C and ^''O 
as a function of the outer mass fraction log(l — Mr/Mt) for our 
0.870 Mq sequence at various selected stages on its cooling brancli, 
as specified by the luminosity and effective temperature values 
(log(L/LQ); logTcff ) in each panel. 

During the PG 1159 and hot white dwarf stages, ^He, 
^^C and ^^O are by far the dominant species in the enve- 
lope of our models. Among the main remaining con- 
stituents are ^^C, and ^^Ne. In particular, as a 
result of the violent H burning during the born-again 
episode, ^'^C reaches abundances as high as 0.05 by mass 
in the envelope. Once He burning becomes virtually 
extinct after the point of maximum effective tempera- 
ture, the H-deficient remnants settle upon their cooling 
track. Here, the chemical abundance distribution will be 
strongly modified by the various diffusion processes act- 
ing during white dwarf evolution. A glimpse of the inner 
chemistry variations that occur during this stage is given 
in Figs. [3] and H for the 0.870 and 0.515 sequences, 
respectively. The upper panel of each figure shows the 
chemical stratification at the start of the cooling track. 
Clearly, gravitational settling causes He to float to the 
surface and heavier elements to sink. Thus, as the evolu- 
tion proceeds, the He abundance in the outer layers in- 
creases, thus turning the PG 1159 star into a white dwarf 
of spectral type DO (Dehner & Kawaler 1995; Gautschy 
& Althaus 2002). At this point, we must mention that 
our treatment of diffusion cannot predict when a PG 
1159 star is transformed into a DO white dwarf. To do 
this, a more elaborated treatment of diffusion in the sur- 
face layers than that we consider here — including mass 
loss through winds, and radiative levitation — should 
be used (Unglaub & Bues 2000). Note the formation of 
a diffusively-evolving double-layered chemical structure, 
that is, a pure He mantle and an intermediate remnant 
shell rich in He, carbon and oxygen result of the last 
He thermal pulse. Because of the more intense gravita- 
tional field and the smaller He content in more massive 
remnants, the double-layered structure in our 0.87 
sequence almost disappears before the DBV instability 
strip is reached at Tcs ~ 30 000 K, see the bottom panel 
of Fig. [31 Indeed, the outer layer chemical profile of the 
massive remnant resembles that of a single-layered struc- 
ture by the time cooling has proceeded to this stage. As 
shown by Fontaine & Brassard (2002) and also by Al- 
thaus & Corsico (2004), the presence of a double-layered 



Fig. 4. — Same as Fig. [3]but for the 0.515 Mq sequence. 

structure bears consequences for the pulsational prop- 
erties of variable DB white dwarfs. In particular, the 
double-layered structure causes the period-spacing dia- 
grams to exhibit mode-trapping substructure (Althaus 
& Corsico 2004). Thus, a complete assessment of the 
evolutionary history of H-deficient white dwarf progen- 
itors like that considered here constitutes a key aspect 
that has to be considered in asteroseismological studies 
of DBV stars. In this sense, we conclude that on the basis 
of a full account of progenitor evolution, a double-layered 
structure is not expected in H-deficient white dwarfs with 
masses higher than about 0.87 Mq by the time evolution 
has proceeded to the domain of the DBV instability strip. 

Our evolutionary sequences are appropriate for mass 
determinations of H-deficient white dwarfs, particularly 
at the early and hot stages when white dwarf evolution 
is still sensitive to the initial conditions. In Fig. [5] we 
show the evolution of our sequences in the log g — logTcs 
diagram and the effective temperatures and surface grav- 
ities (in cm/s^) derived by Hiigelmeyer et al. (2006), 
Dreizler & Werner (1996) and Dreizler et al. (1997) of 
the observed DO white dwarfs. They are listed in Table 
[3] together with the stellar masses determined by these 
authors using the evolutionary tracks of Wood (1995) for 
carbon-core white dwarfs and a He layer mass of 10"''' 
(column 4). In passing, we note that the stellar mass of 
J204158.98-)-000325.4 listed in Hiigelmeyer et al. (2006) 
is not correct, it should read 0.49 instead of 0.6 Mq. We 
have corrected this in Table [31 From Fig. [Sj it is clear 
the importance of finite-temperature effects during the 
early evolutionary stages, particularly in the case of low- 
mass sequences. We expect then that the stellar mass 
for hot DOs with low gravities determined from cooling 
tracks that neglect the progenitor evolution will be dif- 
ferent when our new evolutionary tracks are employed. 
Thus, we redetermined the stellar mass of the DO white 
dwarfs listed in Table [3] using our new sequences. The 
results are given in column 5. For most DOs, mass deter- 
minations from both sets of tracks agree reasonably well, 
but in the case of the hot and lowest-gravity DOs, our 
evolutionary tracks lead to an increase in the derived stel- 
lar mass of up to 10% when compared with the determi- 
nations based on the cooling tracks of Wood (1995). We 
must caution however that some simplifications done in 
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Fig. 5. — Surface gravity g (in cm/s^) as a function of T^fj for our H-deficient sequences with tliick He envelopes. From bottom to top, 
curves correspond to sequences witii stellar masses from 0.870 to 0.515 Mq. We also plot isochrones ranging from 0.03 to 5 Myr measured 
from the highest effective temperature point. The observed DO white dwarfs analyzed by Hiigelmeyer et al. (2006), Dreizler & Werner 
(1996), and Dreizler et al. (1997), as listed in Table [S] are also shown (filled circles and triangles, and star symbols, respectively). The 
thick solid line corresponds to the wind limit for PG 1159 stars taken from Unglaub & Hues (2000), and the two thick dashed lines to the 
wind limits for PG 1159 stars with H/He=0.1 and 0.01 (upper and lower line, respectively). 



the calculations of Wood (1995), in particular regarding 
the chemical profile and He envelope mass of the models, 
prevents us from a consistent and meaningful comparison 
between both sets of tracks. For instance, the decrease in 
surface gravity resulting from the assumption of a pure 
carbon core in the models of Wood (1995) is in part com- 
pensated for by the lower He envelope mass that charac- 
terize these models, as compared to our carbon-oxygen 
core white dwarf models that are characterized by more 
massive He envelopes. Finally, note that for the low- 
gravity DO, J075606.36 -I- 421613.0, we derive a stellar 
mass of 0.467 M0. Although this value results from ex- 
trapolation of our tracks, it is interesting to note that the 
resulting mass is substantially higher than the 0.42 Mq 
quoted by Hiigelmeyer et al. (2006). The new stellar 
mass seems to suggest that the progenitor of this star 
could probably be a He-SdO star formed in a delayed He 
core flash (Miller Bertolami et al. 2008). 

Our sequences provide also realistic evolutionary ages, 
particularly for the hot stages of white dwarf evolu- 
tion, for which employing artificial procedures to gen- 
erate starting white dwarf configurations may lead to 
evolutionary sequences with erroneous thermomechani- 
cal structures and, thus, to wrong ages. We employed our 
tracks to derive evolutionary ages for the white dwarfs 



plotted in Fig. [5l They are listed in column 7 of Table [31 
The zero-point age corresponds to the highest effective 
temperature for each track, see Fig. [2j For the stellar 
masses not covered by our tracks, ages were obtained 
from extrapolation, so they should be taken with some 
caution. We also plot isochrones in Fig. [5] for better 
visualization. Note that the bulk of observed DO white 
dwarfs have ages between 0.6 and 2.5 Myr, though the 
hottest members of the sample are indeed very young, 
with ages less than 0.10 Myr. Note also that the mass 
distribution of young DO white dwarfs differs consider- 
ably from that of the older ones. In fact, for ages less 
than about ~ 0.4 Myr, DO white dwarfs less massive 
than « 0.6 Mq are 75% of the total number of DOs at 
those ages, whilst for longer ages, this contribution de- 
clines to 40%. Finally, the population of high-gravity 
DOs increases considerably for ages greater than 0.4 Myr, 
below Toff « 85 000 K. These features may help to shed 
light on the formation and evolutionary status of DOs. 

In Fig. [7] we show the age (Myr) as a function of 
the luminosity (in solar units) for all the H-deficient se- 
quences. At very high luminosity stages, residual He 
shell-burning constitutes an appreciable fraction of the 
surface luminosity. For instance, at logTcfi = 5, the He 
shell burning luminosity represents 25% of the surface lu- 
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TABLE 3 

Stellar parameters of the DO white dwarf samples of Hugelmeyer et al. (2006), Dreizler & Werner (1996), and Dreizler 

ET al. (1997) — ROMAN AND GREEK LETTERS, AND NUMBERS, RESPECTIVELY. COLUMN 4 GIVES THE STELLAR MASS IN SOLAR UNITS AS 
DERIVED BY THESE AUTHORS USING THE COOLING TRACKS OF WOOD (1995). COLUMN 5 AND 6 LIST, RESPECTIVELY, THE STELLAR MASS IN 
SOLAR UNITS DERIVED FROM OUR EVOLUTIONARY TRACKS WITH THICK AND THIN He ENVELOPES. THE NEXT TWO COLUMNS GIVE THE WHITE 

DWARF AGES. THE EFFECTIVE TEMPERATURE AND GRAVITY OF KPD0005 + 5106 HAVE BEEN TAKEN FROM WERNER ET AL. (2008b). FOR 
THIS STAR, THE LOWER LIMIT TO LOG g HAS BEEN ADOPTED. THE ZERO POINT FOR THE AGES IS AT THE HIGHEST EFFECTIVE TEMPERATURE. 
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minosity in our lowest-mass sequence. This contribution 
declines steeply to 10% at log Toff — 4.95. Except for 
these very high luminosity stages, evolution during the 
early stages is driven by neutrino emission and the release 
of gravothermal energy. In particular, neutrino losses 
exceed photon luminosity during most of the hot white 
dwarf stages. Indeed, except for the hot and low-gravity 
DO sample shown in Fig. [5l it is during this "neutrino 
epoch" that most DO white dwarfs are observed. Note 
from Fig. [7] that the imprints of neutrino emission on the 
cooling curve are more noticeable as the mass increases. 
Less massive sequences are older than more massive se- 
quences for log(L/L0) > — 1. At this stage, neutrino 
emission arrives at its end, thus causing a change in the 



slope in the cooling curve around log(L/L0) ~ —1, and 
later for the less massive sequences. At lower luminosi- 
ties, Coulomb interactions become more important, in- 
creasing the specific heat, until crystallization and the 
associated release of latent heat start at the center (Van 
Horn 1968; Winget et al. 2009). At this point, cooling 
can be well understood on the basis of the simple cool- 
ing theory of Mestel (1952), which predicts less massive 
white dwarfs to cool faster. 

For all our cooling sequences, we computed accurate 
colors and magnitudes based on non-gray LTE model 
atmospheres (Rohrmann et al. 2002). The calculations 
were done for a pure helium composition and for the HST 
ACS filters (Vega-mag system) and UBVRI photome- 
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Fig. 6. — Mass-radius relations for our H-deficient sequences with 
thick He envelopes, corresponding to several effective temperatures. 
From bottom to top we show relations for 40 000 K to 124 000 
K in temperature steps of 4 000 K. In the interest of clarity, the 
mass-radius relations for effective temperatures 40 000 K, 60 000 
K, 80 000 K, 100 000 K and 120 000 K are displayed with dashed 
lines. 




Fig. 7. — Age (in Myr) versus luminosity (in solar units) for our 
H-deficient white dwarf sequences. 

try. Briefly, each model atmosphere is approximated, lo- 
caUy, by plane-parallel layers in hydrostatic equilibrium 
and local thermodynamic equilibrium. Convective en- 
ergy transport in cooler atmospheres has been calculated 
with the so-called ML2 parameterization of the mixing- 
length theory. The correct total energy flux transported 
by radiation and convection is obtained from a Rybicki 
scheme, see details in Rohrmann et al. (2002). Chemical 
populations are computed within the occupation proba- 
bility formalism of Hummer & Mihalas (1988). For he- 
lium pure models, the opacity sources include bound- 
free (He, He+, He2+) and free-free (He, He+,Hc2+, He") 
transitions, electronic and Rayleigh (He) scattering, and 
the most signiflcant line series of He and He+ . The level 
occupation probabilities are explicitly included in the 
evaluation of the line and continuum opacities. As an ex- 
ample of the present atmosphere calculations, we show in 
Fig. [5] the F814W, F606W-F814W color-magnitude di- 
agram for all our sequences. Colors have been computed 
for effective temperatures lower than 40000K, because 
NLTE effects become important above this temperature. 
For higher effective temperatures, most of radiation flux 
emerges at shorter wavelengths. We then expect colors 



to present a markedly less dependence with temperature. 
Comparison between our color sequences and those due 
to Bergeron (see Fig. |8]) reveals a rather good agreement, 
with some small discrepancies presumably due to differ- 
ences of the constitutive physics used in the atmosphere 
codes. 

3.2. Thin envelope sequences 

In section 3.1 we studied sequences characterized by 
the thickest He-rich envelope that can be left after a 
born-again episode. However, there is theoretical and 
observational evidence that He-rich envelopes could be 
substantially thinner. For instance, using asteroseismol- 
ogy, Kawaler & Bradley (1994) derived a He-rich enve- 
lope of ~ 1O~''M0 in at least one of the PG 1159 stars 
with masses 0.6 Mq. More recently, Althaus et al. (2008) 
showed that PG 1159 star models with thin He-rich en- 
velopes appear to be needed to solve the discrepancy 
between the observed (Costa & Kepler 2008) and the 
theoretical rates of period change in the pulsating star 
PG1159— 035. A more direct evidence for the existence of 
PG 1159 stars with thin He envelopes is provided by the 
strong mass-loss events reported in the post-born again 
star V4334 Sgr, also known as Sakurai's object (van Hoof 
et al. 2007). All these facts point at the possibility that 
some PG 1159 stars could be characterized by He-rich 
envelopes thinner than expected from the standard the- 
ory of stellar evolution (Althaus et al. 2009a). Last 
but not least, the mere existence of the hot DQ class of 
white dwarfs (Dufour et al. 2007, 2008) could be reflect- 
ing the existence of He-rich white dwarfs with extremely 
thin He envelopes (Dufour et al. 2008; Althaus et al. 
2009b). In view of these considerations, we examined 
the evolution of H-deficient white dwarfs for the extreme 
situation in which a very thin He envelope could result 
from prior evolution. To this end, we computed five se- 
quences with stellar masses of 0.53, 0.565, 0.609, 0.664 
and 0.741 Mq characterized by He-rich envelopes with 
masses ^ 10~^M*, a value which is of the order required 
by the diffusive/convective mixing scenario to explain the 
origin of hot DQ white dwarfs (Althaus ct al. 2009b). To 
construct initial models with such thin He-rich envelopes, 
wc artificially reduced the He content of our models to 
~ 10~*M* well before they reach the PG 1159 domain. 
These sequences were evolved down to Toff « 20 000 K, 
where the He-rich envelope is diluted by the outer con- 
vection zone. 

In Fig. [5] the resulting tracks in the log 5 — logTcff 
plane are compared with the sequences corresponding to 
thick He envelope models analyzed in section 3.1. When 
thin He-rich envelopes are considered, the evolutionary 
tracks in this plane are strongly modified for the PG 1159 
regime, with important consequences for PG 1159 spec- 
troscopic mass determinations (Althaus et al 2008). For 
the white dwarf regime, models with thin He envelopes 
are more compact — see Fig. [TU] — and, thus, are char- 
acterized by higher gravities. As a result, lower white 
dwarf masses would be inferred when the sequences with 
thin He envelopes are used, as it is clear by examining 
column 6 in Table [31 From Fig. [Slwe expect this to be 
relevant only for the hot white dwarf stages. Indeed, note 
from Table [3] that only for hot DO white dwarfs we find a 
noticeable, albeit small, change in the stellar mass when 
the models with thin He-rich envelopes are used. In Fig. 
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Fig. 8.— F814W versus F606W-F814W color-magnitude dia- 
gram for our DB white dwarf sequences with pure He envelopes. 
Dotted lines correspond to Bergeron et al. (1995) sequences for 
He-atmospheres white dwarfs and stellar masses of, from top to 
bottom, 0.5, 0.6, 0.7, 0.8, and 0.9 Mq. 
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Fig. 9. — Surface gravity as a function of T^g for the sequences 
with thin and thick envelopes and stellar masses of 0.530, 0.565, 
0.609, 0.664 and 0.741 M©. 

111! we show the white dwarf tracks for thin He-rich en- 
velopes together with the observed DO white dwarfs as 
hsted in Table O Note that the low-mass sequences with 
thin He envelopes evolve faster during the hot stages of 
evolution. In fact, for these DO white dwarfs significantly 
smaller ages are inferred when the tracks with thin He 
envelopes are used, see column 8 in Table [31 This is be- 
cause the He shell burning does not represent an energy 
source in the sequences with thin He envelopes. Finally, 
the formation of a double-layered structure in pulsating 
DB stars is not expected if H-deficient white dwarfs are 
formed with sufficiently thin He-rich envelopes. 

4. DISCUSSION AND CONCLUSIONS 

Motivated by the increasing number of detected hot H- 
deficient white dwarfs (of the DO type) as well as by the 
lack of realistic sets of evolutionary tracks appropriate for 
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Fig. 10. — Mass-radius relations for H-deficient sequences with 
thick (solid linos) and thin (dotted line) He envelopes. Curves 
correspond to effective temperatures of, from bottom to top, 40 000 
K to 124 000 K in temperature steps of 4 000 K. 

their interpretation, we computed new evolutionary se- 
quences and colors for non-DA white dwarfs with stellar 
masses of 0.515, 0.530, 0.542, 0.565, 0.584, 0.609, 0.664, 
0.741 and 0.870 Mq. In contrast to previous studies, 
these sequences were obtained considering the full evolu- 
tion of the progenitor stars, from the ZAMS through the 
thermally pulsing AGE phase and born-again episode, 
to the domain of the PG 1159 stars — see Miller Berto- 
lami & Althaus (2006) for details. Our H-deficient white 
dwarf evolutionary tracks are thus not affected by incon- 
sistencies arising from artificial procedures to generate 
the starting white dwarf configurations, and so, they are 
appropriate for realistic mass and age determinations of 
H-deficient white dwarfs, particularly for objects at the 
early and hot stages. Calculations were followed to ad- 
vanced phases of evolution down to an effective temper- 
ature of 7000 K. 

We applied this new set of tracks to redetermine 
the stellar masses of all known DO white dwarfs with 
spectroscopically-determined effective temperatures and 
gravities. The results are listed in Table \3\ The stel- 
lar masses for hot DOs with the lowest gravities inferred 
from our tracks differ appreciably from the determina- 
tions based on the evolutionary tracks of Wood (1995). 
This is because these tracks do not consider the full evo- 
lution of progenitor stars. This is also true for the age 
determinations of these stars. In particular, our evo- 
lutionary sequences provide realistic ages at very early 
stages of evolution. For thick He layers, we find that 
He shell-burning supplies an appreciable fraction of the 
surface luminosity, thus changing appreciably the evolu- 
tionary time scales, particularly of the hottest DO white 
dwarfs with low gravities. The amount of He left in the 
remnant white dwarf affects both the mass and the age 
determinations of hot DOs. Models with thin He-rich en- 
velopes are characterized by higher gravities than models 
with thick He-rich envelopes, thus predicting lower white 
dwarf stellar masses for the observed DOs. 

The evolutionary tracks presented here cover all the 
stages from the domain of luminous PG 1159 star region 
to the hot white dwarf domain. Thus, they constitute 
an homogeneous set of evolutionary tracks to study both 
type of stars. This is of major interest because we can 
make consistent mass determinations for both the DOs 




Fig. 11. — Same as Fig. [5]for our H-deficient sequences with thin He envelopes. We also plot isochrones (Myr) ranging from 0.03 to 5 
Myr. The location of the observed DO white dwarfs listed in Table[3]is also shown. 



and PG 1159 stars for the first time. In Fig. [T^lwe com- 
pare the mass distribution of the DO white dwarfs hsted 
in Table [3] with the mass distribution of PG 1159 stars 
constructed on the basis of mass determinations made by 
Miller Bertolami & Althaus (2006), shaded region and 
solid line, respectively. The size of the bins in both dis- 
tributions is 0.1 Mq. Note that massive DOs outnumber 
massive PG 1159 stars. We obtain a mean DO mass of 
0.644 Mq, considerably higher (by 0.071 M©) than the 
mean mass of PG 1159 stars, 0.573 Mq. One may wonder 
if the reason for this difference is simply that evolution 
through the PG 1159 stage proceeds considerably faster 
for massive stars (Miller Bertolami & Althaus 2006), thus 
decreasing the possibility of detection. To elaborate on 
this point, we have proceeded as follows. In order to 
avoid the possibility that the massive PG1159 stars are 
under-represented because they do not spend long times 
in the PG 1159 stage, we have tested the hypothesis that 
PG 1159 stars and DO white dwarfs form an isolated 
evolutionary sequence by analyzing the PG 1159-l-DO 
group as a whole. Taking advantage of the consistent 
set of tracks for both groups, we have divided the sam- 
ple (78 stars) into two groups based on their age. The 
young group, which consists of all stars with ages below 
0.3 Myr (this includes all 37 PG 1159 stars plus 12 young 
DO white dwarfs) and the old group consisting of all DO 
white dwarfs with ages greater than 0.4 Myr (29 stars). 
It must be noted that there are no PG 1159 stars in the 



old group. 

The choice of the 0.3-0.4 Myr isochrone to divide both 
groups is admittedly arbitrary but suggested by the lo- 
cation of the different wind limits studied by Unglaub & 
Bucs (2000). Dividing the whole PG 1159+DO group by 
an isochrone allows to test the hypothesis that PG 1159 
and DO stars form an isolated evolutionary sequence. 
This can be done comparing their mass distributions (as 
relative numbers within a given group will not be af- 
fected by differences in evolutionary speeds). The two 
groups show significant differences in their mass distri- 
butions. While the young group has a mean mass of 
My = 0.584 Mq (with a variance of cry = 0.080 Mq) the 
old group has a mean mass of Mq = 0.655 Mq (with a 
variance of ao — 0.087 Mq). Thus, the difference in the 
mean masses is comparable with the variance within each 
sample, suggesting that both mass distributions are sig- 
nificantly different. This is confirmed by a Kolmogorov- 
Smirnov test. The result of this test is that the proba- 
bility that both samples are taken from the same mass 
distribution is rather small (~ 10~^). 

These results clearly suggest that there may be some 
other important evolutionary channels operating within 
PG 1159 and DO stars. For instance, it is possible that 
not all PG 1159 stars evolve to DOs. In particular, those 
PG 1159 stars resulting from a LTP — low-mass rem- 
nants are more prone to experience a LTP episode — 
are expected to evolve into DA white dwarfs (near the 
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Fig. 12. — The shaded histogram shows the mass distribution 
for the DO white dwarfs Usted in Table 3. The mass determina- 
tions have been done using the evolutionary tracks described in this 
work. The solid line corresponds to the PG 1159 mass distribution 
from Miller Bertolami & Althaus (2006), and the dashed line to 
that of DB white dwarfs taken from Kepler et al. (2007). 

winds limits shown in Fig. [5] with sohd dashed hnes), 
thus avoiding the DO stage, as shown by Unglaub & 
Bues (2000). This is so because, in contrast to a VLTP 
event, in a LTP, H is not burned, but instead diluted 
to low surface abundances (Miller Bertolami & Althaus 
2006). Alternatively, some DO white dwarfs could result 
from evolutionary channels that not involve only the PG 
1159 stars. For instance, they could be the result of post- 
merger evolution involving the giant, H-deficient RCrB 
stars, via the evolutionary link RCrB ^ EHe (extreme 
He stars) He-SdO+ ^O(He) ^ DO — see Ranch et 
al. (2006) for a connection between 0(He) stars and 
RCrB stars. This is reinforced by the recent study by 
Werner et al. (2008a,b) of the star KPD0005-I-5106, the 
hottest known DO white dwarf. These authors present 
evidence that KPD0005-I-5106 is not a descendant of PG 
1159 stars, but more probably related to the 0(He) stars 
and RCrB stars. 

The group of young and hot DO white dwarfs is partic- 
ularly relevant. Note that they are located considerably 
above the wind limit for PG 1159 stars (see Fig. [5]). In 
fact, above this line, mass-loss is large enough to prevent 
gravitational settling. Hence, the transformation of a 
PG 1159 star into a DO white dwarf should be expected 
approximately below this line. Although the estimated 
mass-loss rates could be overestimated by more than a 
factor often — in which case the transition from PG 1159 
to DO would be near a line with \ogg w 7.5 (Unglaub & 
Bues 2000) — it could be possible that the less massive 
and young DOs are not the descendants of PG 1159 stars. 
This is also suggested by the fact that for ages less than 
about 0.3-0.4 Myr, the DO white dwarf population is 
markedly less massive than at greater ages. Indeed, the 
number of detected old, less massive DOs is smaller than 
the number of young, less massive DOs, despite the evo- 
lutionary timescales being a factor about 5 for the former. 



This could be indicative that some of the less massive 
DOs experience a tranformation in their spectral type 
after w 0.3 — 0.4 Myr of evolution. In particular, traces 
of H could turn the white dwarf spectral type from DO to 
DA as a result of gravitational settling. An approximate 
location where this transformation should occur is pro- 
vided by the thick dashed lines shown in Fig. [5l which 
provides the wind limits for PG 1159 stars with differ- 
ent initial H abundances (Unglaub & Bues 2000). It is 
worth mentioning that because of the much smaller car- 
bon content and correspondingly less intense winds, such 
limits would be expected at lower gravities in a DO white 
dwarf with similar H content. According to Unglaub & 
Bues (1998), for a DO white dwarf with logg = 7 and 
H number ratio of 0.01, the transformation into a DAO 
in presence of a moderate mass loss would take about 
0.25 Myr, not very different from the age range where 
the DO mass distribution seems to change. At this point 
it is worthwhile to comment on the possibility that the 
origin of low-mass, hot DOs can be connected with the 
merger of two low-mass. He-core white dwarfs (Guerrero 
et al. 2004; Loren-Aguilar et al. 2009), which will evolve 
to become a low-gravity extreme He star and then a hot 
subdwarf (Saio & Jeffery 2000). The presence of H ex- 
pected in this case, could lead to a transformation of the 
DO white dwarfs into DA ones. 

In Fig. [12] we have also included with dashed line the 
mass distribution of the DB white dwarfs in SDSS hot- 
ter than Toff = 16 000 K taken from Kepler et al. (2007). 
The DB distribution is more weigthed to massive stars. 
In fact, the DB mean mass is of 0.711 Mq, consider- 
ably higher than the mean DO mass of 0.644 M©. This 
marked difference in the mass distribution between both 
populations is certainly quite unexpected in view of our 
understanding that DO stars must evolve into DB stars. 
However, it is possible that the DB stellar masses consid- 
ered in Fig. [T2] could be overestimated. This hypothesis 
was raised by Kepler et al. (2007), who discuss that an 
increase in the measured surface gravity of DBs below 
Toff ~ 18 000 K could result from missing or incorrect 
physics in model atmospheres at such temperatures. It 
could be also possible that the increase in the gravity 
of the DBs by T^e « 18 000 K could be reflecting the 
existence of DBs with traces of carbon in their atmo- 
spheres. As shown by Dufour et al. (2005) in the context 
of cool DQs, the surface gravities derived from model at- 
mospheres including carbon are significantly lower than 
the gravities obtained from pure He models. The idea 
that some carbon could be present in the atmosphere of 
some hot DBs is particularly attractive if a large fraction 
of DB white dwarfs evolve through the hot DQ stage. In- 
deed, as discussed in Althaus et al (2009b) in the diffu- 
sive/convective mixing picture for the hot DQ formation, 
shortly after the formation of the carbon-rich convective 
envelope, by Toff « 20 000 K, a substantial depletion of 
carbon is expected in the entire convective envelope, be- 
cause of the very short diffusion timescale at the base of 
the convection zone. Hence, it is expected that the hot 
DQ stage, during which the white dwarf is characterized 
by high amounts of carbon in its surface, is indeed a 
short-lived phase. However, the change in composition 
in the convective envelope will lead to changes in the 
depth of the convection zone, which will affect, in turn, 
the timescales of diffusion. The lack of detailed computa- 
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tions of this feedback between convection and diffusion 
shortly after the formation of a DQ does not allow to 
make reasonable predictions, but it is not discarded that 
the star experiences additional mixing episodes, eventu- 
ally reaching some stationary situation in which the star 
resembles a He-rich DB white dwarf with some carbon 
in its atmosphere. This being the true course of events, 
it would be worthwhile to re-determine the mass distri- 
bution of DBs in the frame of model atmospheres that 
incorporate the effects of the presence of carbon. Fi- 
nally, in view of the material discussed in this section, it 
is feasible that a fraction of low-mass DOs experiences a 
change in their spectral type at young ages. 

In closing, because our DO tracks consider the evolu- 
tionary history of progenitor, we can make sound pre- 
dictions about the chemical profile expected in the enve- 
lope of H-deficient white dwarfs. This is of relevance for 
asteroseisniological inferences about pulsating DB white 
dwarfs. In particular, we find that depending on the 



stellar mass, diffusion processes lead to a variety of en- 
velope chemical profiles, from single-layered to double- 
layered structures, by the time evolution has proceeded 
to the domain of the variable DBs. Detailed tabulations 
of our evolutionary sequences are available at our web 
site http : / / www . f caglp . unlp . edu . ar/ evolgroup 
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